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1 Introduction 



Future e^e^ colliders will likely have the option of operating in 77 or 67 collision modes [0]. 
These modes are reached by Compton scattering laser light off one or more of the incoming 
fermion beams, and then colliding the resulting high energy photons with the remaining 
fermion beam or with each other. There is a large potential for 67 and 77 collisions to elu- 
cidate possible physics beyond the Standard Model; previous investigations have focused on 
anomalous couplings 0, searches for extra dimensions 0, properties of supersymmetry [^], 
and a broad host of other topics. 

One subject that deserves further study is CP-violating gauge boson self- couplings. 
CP violation is one of the most poorly understood aspects of the Standard Model (SM). 
Present data merely fixes the value of the CP violating phase in the CKM matrix, and cannot 
test if this phase constitutes the only source for CP violation. In fact, studies of baryogenesis 
within the SM suggest that additional CP violating terms are required in order to generate 
the observed baryon asymmetry in the universe |^. Most discussions of CP violation at 
photon colliders work within the context of models such as supersymmetry, and focus upon 
either Higgs or top quark production [Q. However, since almost every extension of the SM 
contains new CP violating phases, it is desirable to eliminate any possible theoretical bias 
by studying CP violation within the generic context of an effective Lagrangian approach, 
without assuming any underlying mechanisms. A handful of such works have been performed 
in the past |[, concentrating on 77 ^ i7, 77 — W'^W~ , or top quark production. Here, 
we extend these studies by examining possible CP violating quartic gauge boson couplings. 
77 colliders are particularly suited to studying such couplings; they can be probed in 2 — >^ 2 
scattering processes, unlike at e+e^ colliders. 

The SM contributions to the processes 77 — >■ 77, 77 'jZ, and 77 — ^ ZZ, including 
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electroweak contributions, were first computed in B. Recently, they have been reexam- 



ined |TT], |T2|, |T3[, and shown to exhibit several interesting features that motivate our 
analysis. The SM amplitudes vanish at tree level, and may therefore be quite sensitive to 
the effects of new physics. At one loop they acquire large imaginary parts, completely dom- 
inated at high energies by the helicity amplitudes A4±±±±, ■M±^±^, and ^A±zpzp±, where 
± denotes the direction of transverse polarization relative to the beam direction. Such am- 
plitudes can interfere strongly with CP violating contributions from new sources, yielding 
greatly enhanced sensitivity to these new effects. In addition, the ability to polarize both 
the initial laser light and fermion beams allows the construction of observables which are 
sensitive only to these interference effects. These properties, together with the experimental 
cleanliness of 77 77, •yZ, and ZZ scattering, suggest that these processes might provide 
powerful tools in searches for new physics. 

In this paper we will show that 77 — > 77, 7Z, and ZZ at a photon collider can 
provide sensitive tests of CP violation in the gauge boson sector. We limit our study to 
genuinely quartic gauge boson operators, as contributions to the two point functions are 
strongly excluded by current data and the three point functions are likely better tested 
elsewhere. We consider manifestly SU(2) x U(l) invariant operators constructed from the 
appropriate field strength tensors, making no assumptions as to the origin of electroweak 
symmetry breaking. As a result, our operators are of dimension eight. Throughout the paper 
we will identify features of our analysis that might be useful for other new physics searches 
using these processes. 

The paper is organized as follows. We first present the density matrix formalism for 
photon-photon collisions in some detail, both for completeness and to motivate the expression 
for the ensemble average of the scattering amplitude, which will be used throughout our 
analysis. We then briefly discuss the construction of the anomalous CP violating operators 
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we have studied. Our results for the three processes considered are presented next, including 
discussions of the various initial laser and fermion beam polarizations. Finally, we present 
our conclusions, including a comparison of our resulting sensitivity to these effects with 
others in the literature. 

2 77 Collisions 

In this section we present the density matrix formalism for photon-photon collisions. Our 
discussion is similar to that found in the literature (see, for example, 0), but is included 
here to provide motivation for eq. (9), which is used throughout our analysis. 

Consider a photon moving along the z axis; its polarization vectors for positive and 
negative helicities are given by 

e^ = ^(0,Tl,-^,0). (1) 

We will denote the photons corresponding to states with these polarizations by |±). The 
most general photon state can be written in this basis as 

= -cos(a)e^"^|+) + sin(a)e-^"^|-) , (2) 

where < a < 7r/2 and — vr < < tt. Writing this state in 4- vector form demonstrates that 
the choices a = and a = 7i/2 lead to circularly polarized states, while a = n/A leads to 
a linearly polarized state. The angle cj) describes the direction of linear polarization in the 
X — y plane. 

The experimental setup in 77 colliders has been described in detail elsewhere [0; 
briefly, one Compton scatters laser photons off fermion beams, and then collides the backscat- 
tered photons. When observing the scattering process 77 X, the energies or helicities 
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of the incoming photons are not known; they will be statistically distributed, with the spe- 
cific distribution being determined by the details of the Compton scattering process. When 
measuring the differential cross section, one no longer determines | {(p finai\-M.\(f)initiai) but 
rather the ensemble average of M.: 

\M\l^,^j:w'\{<l>f,nal\M\<PU,,J\\ (3) 
i 

where the sum over i sums over all possible initial states weighted by w\ the probability of 
their occurrence. We can write this in the form 

I M |L= J2ihinal\M\\j){\k\M\<p final) Pjk , (4) 

where we have introduced the density matrix for the incoming photons in the basis defined 
by Xf 

p,, = ^^^(A,|C)(C|A/c) . (5) 

i 

Let us compute the density matrix for the laser photon of eq. (2); introducing the notation 
Pc — cos(2q;) and Pt — sin(2Q;), we have 

If 1 + Pc -Pte^^^ 



P-2\-Pte-^''^ l-PcJ- 

Pc and Pt measure the amounts of circular and linear polarization, respectively. This is 
consistent with the polarizations for the values a — 77/2 and a — tt/4 noted before. After 
the Compton scattering process, the degrees of linear and circular polarization are described 
by distributions which are dependent upon the fraction of the fermion beam energy the 
laser photons acquire. Denoting this fraction by x, and the helicity distribution functions 
by Cc,t{^)j density matrix becomes 



The eigenvalues of this matrix are 



(8) 



we no longer have a pure state unless + ~ ^■ 

As the processes considered in this paper involve two Compton scattered photons, the 
complete density matrix will be the tensor product of the density matrix for each photon. 
Also, since one of the initial laser photons will be moving along the —z axis, we must take 
(f) —(f) in its density matrix. Referring back to eqs. (4) and (7), we can write the ensemble 
average for the process 77 — > X as 



M 



^ I I M++ I' + I M— I' + I M+- I' + I M- 



+ec(a;i) I M++ I' - I M— I' + I M+- I' - I M-+ 



+^0(^2) I M++ I' - I M— I' - I M+- I' + I M-+ 



+ic{xi)ic{x2) 



M++ r + 1 M— r - 1 M+- \-\M- 



■+ 



I* c2«</'l 



-2^t{xi) Re [M++M%e''^' + M+-M*__e 



+2^t{xi)Ux2) Re M+-M*__e'''l'' - M++M*e^'^' 



+2Uxi)Ux2) Re \M—M%e''">>' - M+-M*e^'^' 



(9) 



where, for example, denotes the helicity amplitude for incoming photons with helicities 

+1 and —1. We have suppressed the final state X; depending upon the process considered 
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the observable final state might require a sum over various helicities. This expression will 
be used to construct observables for all the processes considered in this paper. 

The physical cross section involves a convolution of this scattering amplitude with 
the photon number density function of each photon. The explicit form of this function, as 
well as the forms for the linear and circular helicity distribution functions, can be found in 
the Appendix. 

3 Constructing CP Violating Operators 

Here we construct the most general set of operators that contribute to neutral gauge bo- 
son self-interactions, subject to the following constraints. We consider only SU(2) x U(l) 
invariant operators, and further restrict these to CP-odd terms. As the effects of CP-odd 
trilinear operators have been extensivley studied we consider only those terms that 



lead to quartic or higher gauge boson interactions. We also make no assumption as to the 
mechanism of electroweak symmetry breaking. These restrictions lead us to the following 
set of seven dimension eight operators constructed from the W and B field strength tensors: 

O^BB)(BB) = {B^,B^^){BP'^ Bp.). 0(^WW)(WW) = {WapuWniWr Wbpa), 
OiBB)(WW) = iB,,,Bn{Wr Wap.), 0^wW)iBB) = {Wa^^WDiB'^ Bp.), 
O^BW)iBW) = {B^^WDiB''' Wapa), OwBWB = B^.W^ B.p, 

ObWBW = B^^'WaupB^^ Waa^. • (10) 

Here W^'^ is the SU(2) field strength tensor and P^*^ the U(l) field strength tensor. We have 
introduced the notation 

~ 1 

Bp,u= - ^pupa B'^^ , (11) 



and our convention for the e tensor is 60123 = 1- A brief comment on the completeness of 
these operators is in order. CP violating operators containing three e tensors are reducible 
to those containing only one such tensor, as the product of two e tensors can be written as 
a determinant of metric tensors. Operators such as e^^i^p^A'^'^B'^C'^^D'^, where the e tensor 
contracts an index on each field strength tensor, are reducible to those where the e tensor 
fully contracts one of the field strength tensors through use of the Schouten identity (see the 



first reference in [jT^, |jT5| 



+ 9ape 

which states that no tensor antisymmetric in five indices exists in four dimensions. Similarly, 
the operators 

WrWaupBP'' 5.,., Bf^'B^pWr Waap , (13) 

can be reduced to those listed in eq. (10) through use of the identity 

epo.p.F^'^F^'^ = ^gi^e^p.sF'^f'F^' . (14) 

This is true for F^^" antisymmetric in four- dimensional spaces (an analog holds in arbitrary 
dimensions [|16[)- Finally, since the processes considered here only involve W^'^ , the operators 

{WapuWDiWr Wbpa), WrWt^pWr Wba^. (15) 

become equivalent to the second operator in eq. (10) through use of the above identity, while 
the operators 

eabciWap.WniB''' Wcpa), e^UB p.WniW^ Wcpa), 

eabcWrWku.Wr B.p, e.bcB^'-Wa.pWr Wcap (16) 



only contribute to processes involving W^. We thus see that the seven operators delineated 
in eq. (10) form a complete and independent set. 

Each of these operators will give rise to a number of different 7777, 777 -Z^, and 77ZZ 
operators, as well as other quartic operators involving three or more Z bosons and various 
quintic terms. In this paper we will concentrate on those quartic operators relevant for the 
scattering processes 77 — > 77, 7Z, and ZZ; the other terms cannot be probed in 2 — > 2 
scattering, and the resulting constraints obtainable on the operators of eq. (10) would be 

weakened. Denoting the operators to be studied by O]'^ , O]^ , and Of^, each SU(2) x U(l) 
operator will have an expansion of the form 

osmxud) ^ ^i^or + h\of + ejD^ + . . . , (17) 

where the a, 6, and c are functions of the weak mixing angle and the ellipsis denotes the ne- 
glected quartic and quintic terms. To determine the independent 77, 7Z, and ZZ structures, 
we set 

W^"" = sin(%) F'^'^ + cos(M , (18) 

where 

Ffiv — d^iAy — di^Afj^ , 

= di,Z,-d,Z^. (19) 

Doing so, and making use of the identity in eq. (14), leads to the following forms for the 
Of: 

Of_^ : {F^,F^-){F'^'^Zp.),{F^,Z^-){FP'^Fpa),F^^F,,FP'^Za^, 



F^'^Z^pF"'' , Z^-F^pZ^'^ F.^ . (20) 

We will present the coefficients a, 6, and c when we discuss the relevant process. 

We will organize our study as follows. A section will be devoted to each of the 
processes 77 77, 77 7Z, and 77 ZZ. In each section we will discuss the properties 
of the relevant Of ^ focusing on the observable asymmetries associated with each operator and 
the role of varying the initial polarization states. We will then discuss the measurement of 
the various (9^f^(2)x(7(i)^ g^j-^^ estimate the sensitivity to these operators that can be obtained 
at a photon collider. We will bound each of these operators separately; each is multiplied 
by an arbitrary coefficient that could cause contributions from different operators to cancel, 
but we will ignore this possibility. To preview what we will observe, it will turn out that 
in certain (9^^(2)xc/{i) ^j^g structures of eq. (20) will interfere destructively, while in others 
they will combine constructively, leading to widely varying sensitivity to the 0^f^(2)x'^(i)_ 



4 77 ^ 77 

In this section we consider the process ^{k\) + 7(^:2) liPi) + 7(^2)- We also present and 
motivate the various parameterizations and approximations used throughout the paper. A 



detailed presentation of the Standard Model amplitudes can be found in |I0|, [Tl|]; we focus 
here only upon those features relevant to our analysis. 

Let us ffist discuss the properties of the CP violating amplitudes we are considering. 
The detailed expressions for these amplitudes, and those for the other processes considered 
in this paper, can be found in the Appendix. Bose symmetry implies the relations 

M^,Us, t, u) = MZUs, u, t) = MZcis, t, u) , (21) 



where s, t, and u are the usual Mandelstam invariants: s = {ki + A;2)^, t = {ki — piY, and 
u = {ki — ^2)^- We have denoted the CP violating amplitudes with the superscript CP. The 
operators considered in this paper are odd under both P and T, which implies 

M'^Us, t, u) = -M^^_,_,_,is, t, u) = -A^S:,(s, t, u) . (22) 
In addition, the standard crossing symmetries hold. These relations immediately imply 

MZ±± = = ^±TT± = • (23) 

The only non- vanishing amplitudes are ■M.±±zfz:f The Standard Model amplitudes are domi- 
nated at high energies by the amplitudes M±±±±, -^±^^±=p) and M.±^:f±- These amplitudes 
are primarily imaginary, although the real parts are non- negligible; features of all 77 77 
amplitudes are discussed in [0. All of the CP odd asymmetries that can be constructed 
involve the interference of Ai±±^zfi with a Standard Model amplitude; to increase our sen- 
sitivity to the anomalous interaction we should attempt to define an observable with inter- 
ference between A^^l^^^ and one of the dominant Standard Model terms. Unfortunately, 
no observable exists which contains such an interference with one of the large imaginary 
SM amplitudes. We can, however, construct the following observable which contains an 
interference with one of the real pieces of the SM amplitudes: 



1 1 


?01 d(j)2 


\dn J 


5(01 - 02 


- ^/4) - ) 5(01 - 


02 + vr/4) 


Jo Jo a 


l(t>i d(j)2 


Y dcr^ 

\dn ) 


5(01 - 02 




02 + 7r/4) 



(24) 



where 0i, 02 denote the angles of linear polarization introduced in our discussion of the den- 
sity matrix formalism. The experimental determination of this asymmetry requires several 
measurements, as indicated in eq. (24) by the integration over linear polarization directions. 
The delta functions in A^^ denote the need to maintain a fixed angle between the linear 
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polarization directions of the initial beams. This observable has been considered previously 
in J^, 0. Referring to eq. (9), and implementing the relations in eqs. (21-22), we can show 
that the numerator of this asymmetry contains the term 

A,, oc Uxi)Ux2) Re (-Mf+'++) Im • (25) 

As shown in the Appendix, the CP violating amplitudes are purely imaginary, leading to a 
non- vanishing result. As this observable is sensitive only to the degree of linear polarization, 
we will use the following set of initial parameters: 

Pcl = Pc2 = Q] Pti=Pt2 = i; (26) 

we note that with no circular polarization the fermion beam polarizations do not enter any 
expression. We can simultaneously measure another independent quantity, the denominator 
of Ay^. This is the unpolarized cross section 



M++ |2 + I M— I' + I M+- |2 + I M-+ , (27) 



dcos{6) 2567rs 

the form of which can also be obtained from eq. (9). 

We will now consider the operator relevant to 77 final states in eq. (20). We write it 
in the form 

Or = jf^{F,.Fn{F''^F,.), (28) 

where Ai is the energy scale of the physics leading to these operators, and Qeff the associated 
coupling constant. We will hereafter set Qeff = 1, effectively absorbing this coupling constant 
into the definition of Ai. Although we will quote our sensitivity in terms of the energy scale 
that can be probed, the reader should realize that this is not exactly the scale associated 
with the new physics leading to these anomalous couplings. We have also ignored a possible 
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±1 appearing in front; the unpolarized cross section is unaffected by this sign, and A^^ only 
acquires an overall sign change. 



In our analysis, we have used the approximate SM amplitudes found in W^, valid for 
m^r/{s,t,u} < 1, which holds for the collider energies considered here. We will use similar 
approximate expressions when considering 77 7Z and 77 ZZ. The features of the SM 
amplitudes noted earlier occur when the energy of the process becomes greater than several 
hundred GeV. In accordance with this approximation we employ the cuts 

I cos{e) \< 0.866, VOA < X, < x^ax ■ (29) 

Throughout this paper, we take the e+e~ center of mass energy to be ^/s = 1000 GeV 
when obtaining our search reaches. In this case, at the endpoints of this region, where 
xi = X2 = VOA and cos(6') = 0.866, we have m^/ | t |~ 0.24; this value becomes much 
smaller as we approach the center of the parameter region, thus validating our approximation. 
Although a more detailed analysis taking into account detector properties and experimental 
cuts would need to use the complete SM expressions, our study captures the salient points. 

We show the total cross section and asymmetry versus ^/s for this operator in fig. (1) 
for a "typical" value of Aq, for purposes of demonstration. In presenting these results we have 
assumed an e^e~ integrated luminosity L = 500 fb~^, which is the quoted yearly integrated 
luminosity for planned linear colliders [|l^ . However, a determination of A^^ requires several 



separate measurements, and it is certainly not expected that colliders will operate primarily 
with purely linearly polarized beams, or even primarily in a photon collision mode. These 
numbers should be interpreted as results coming from several years of operation and are 
intended for purposes of illustration. To provide a more conservative outlook, our quoted 
sensitivity will be given as a function of integrated luminosity, beginning at the modest value 
L = 50fb~^ We have also assumed the value a = 1/137.036 in this and all other numerical 
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results presented. 

At the luminosity L — 500 fb~^, the asymmetry is much smaller than the associated 
errors, and only a change in the total counting rate is statistically observable. Such an effect 
can arise from a variety of sources, and the identification of CP violation requires higher 
luminosities to observe a non- vanishing asymmetry. 

The transcription of these results into statements about the SU(2) x U(l) operators 
is simplified by there being only one 77 operator. The coefficients cisjj^2)xu{i) 

"■{BB){BB) — ' "■iWW)(WW) — ' 
^{BB)(WW) — ^{WW)(BB) — (^{BW){BW) — ^W^W i 

1 1 -'-2 2 ('i.f\\ 

^BWBW — ^WBWB — '^^W^W ' W'-'J 

where cw,sw — cos(^vk), sin(6'vy). The SU(2) x U(l) operators of eq. (10) take the form 

Oa^^,{F,^Fn{F'"'F,.). (31) 

To estimate the value of Aq, that can be probed at a 77 coUider we have performed a 
combined least-squares fit to the total cross section and asymmetry. We have assumed 
standard statistical errors and an additional 1 % luminosity error in the integrated cross 
section. The fit was performed with i/s — 1000 GeV; approximate results for other energies 
can be obtained by scaling these numbers. As can be seen from eq. (30) there are only four 
different results for this process; these are presented in fig. (2). 

Although this asymmetry is not particularly sensitive to CP violation as parameter- 
ized by the operators we consider, we note that it is a good test of any CP violating amplitude 
for which any of the real amplitudes A1±±=p=p, -^i^iip; oi' -^i^iFi ^-^^ non- vanishing. In 
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these cases terms of the form 

6(xi)6(x2) Im (^^^^) Re (^^^) (32) 
appear in the fifth hne of eq. (9), leading to a large A^^. 



5 77 ^ "yZ 

In this section we study the effects of anomalous CP violating operators in 7(/ci) + 7(^:2) - 
7(pi) + Z{p2). A detailed study of the SM amplitudes for this process can be found in [|12 

As in 77 — > 77, the jZ amplitudes satisfy certain relations dictated by Bose symme- 
try; the anomalous couplings considered satisfy additional relations because of their trans- 
formations under P. Let Xz = 0, ±1 denote either longitudinal or transverse Z polarizations. 
Bose symmetry implies 

MZxA^,t,u) = i-lY-''Mfj:^^is,u,t), (33) 

and oddness under P requires 

M''Zx,is,t,u) = i-l)^-M^:^_,_,^,^is,t,u) . (34) 

Since the initial and final states are different in this process, the transformation properties 
under T do not imply any relations between the amplitudes. The explicit computation of 
the amplitudes is presented in the Appendix. Unlike in 77 — > 77, there are several surviving 
amplitudes. 

The structure of the SM amplitudes for 77 7Z is similar to that for 77 77, 
and even more pronounced. At high energies the process is dominated by the imaginary 
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Figure 1: Total cross section (top) and asymmetry (bottom) for the anomalous 77 operator, 
with Ai = 2 TeV and L = 500 fb~^. The bars correspond to the statistical error plus a 1 % 
luminosity uncertainty in the event rate; the solid line represents the SM event rate. 
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Figure 2: Sensitivity to for each SU(2) x U(l) operator as a function of integrated 
luminosity at the 95% CL. 

parts of the amphtudes ■M±±±±, ■M±^±zf, and ■M±^:f±', all other amplitudes are completely 
negligible |T2[. As before, we must construct an observable that contains an interference of 
a CP odd amplitude with one of the dominant SM terms. The asymmetry considered in the 
previous section contains only interference with the real parts of the SM amplitudes, and is 
therefore unacceptable. For interference effects between imaginary CP odd and CP even 
amplitudes, we must consider the following asymmetry which is measurable with circularly 
polarized beams: 




This asymmetry has been considered previously in 0. The subscripts ± denote the initial 
polarization states of the laser and fermion beams, which we will now discuss. Setting 
Pti = Pt2 = 0, we are left with four parameters describing the inital state polarization: 







A(BB)(BBjJ 


- 
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— " "'^wm - 




I.J . . . . 1 . , . 
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Pel, Pe2, Pii, and Pi2- We will set \Pe\ = 0.9 and |P/| = 1.0, consistent with the expected 
capabilities of future facilities |T^, and label our inital states as {Pei, Pn, Pe2, Pi2)- In the 

SM, there are six independent states: (+ + ++), (+ + H — ), (+H ), (H 1 — ), ( — hH — ), 

and (H ), where, for example, (H 1 — ) means Pel = 0.9, Pn = —1.0, Pe2 = 0.9, and 

Pi2 = —1.0. States obtained by an overall sign flip are identical; for example, (H \ — ) and 

( — I h) lead to the same observables. This is not true when CP violating interactions are 

present. The asymmetry of eq. (35), where the subscript + refers to a given inital state and 
— to the state obtained by flipping the signs of the polarizations, will vanish in the SM and be 
non-zero in the presence of the anomalous couplings. In addition to this asymmetry we can 
simultaneously measure another independent quantity, the denominator of the asymmetry. 
Referring to the ensemble average in eq. (9) and using the relations in eqs. (33-34), we see 
that this denominator is just twice the polarized cross section. 



da. 



pol 



dcos{e) 1287rs 



{1 + U^i)Ux2)) (\ M 



(36) 



The operators relevant to 7Z final states can be found in eq. (20); we will write them 
as follows: 



of 



of 



Of 



(A 



(P^,P'^")(P'"^Z 



pa) 



^^{F.^ZniF"' Pp.) 



(Aa)^ 



(37) 



Our philosophy in parameterizing these operators is identical to that discussed in the 77 
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77 section; there is again a coupling constant geff associated with each operator that we have 
set equal to 1, and a possible factor of ±1 appearing in each operator that we have ignored. 
The detailed amplitudes are presented in the Appendix; here we will discuss several properties 
relevant to our analysis. The first is that the amplitudes for longitudinally polarized Z states 
are suppressed relative to those with transverse Z states by a factor of mz/y/s, as dictated 
by the Goldstone boson equivalence theorem. This, and the fact that the tagging efficiency 
for longitudinal states is at best 10% , means that we will not explore the effect of selecting 
final state polarizations in our analysis. 

Before presenting our results, let us examine A^z in more detail. Our discussion will 
illustrate the complementary information that can be obtained from examining all of the 
possible initial polarizations. Concentrating on the numerator of Ajz, we see that 



AjZ OC iCcixi) + Cc{x2)) 



M++ \-\M- 



+ iU^i) - U^2)) [I M+- \' - I M-+ \'\ , (38) 

where a sum over final state helicities is implied. The first term is symmetric under the 
interchange of the two initial photons, which implies symmetry under the transformation 
cos(6') ^ —cos{9) and under the interchange of the initial polarization states, {Pgi, Pn) ^ 
{Pe2,Pi2)- The second term is antisymmetric under both of these exchanges. By selecting 
an appropriate initial polarization state, we can isolate each term. For example, the choice 
(+ + ++) is symmetric under the interchange of the two initial polarizations; with this choice 
only the first term of A^z will contribute, whereas the choice ( — h H — ) is antisymmetric 
under interchange of the two initial polarizations, and with this selection only the second 

term contributes. A choice such as (H ) is of mixed symmetry, and is sensitive to both 

terms in A^z- Such a variation of the inital states allows us to isolate the various anomalous 
amplitudes which contribute to A^z is sensitive; with symmetric initial polarizations it is 
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sensitive to A^^^_|__|_, while with antisymmetric polarizations it is sensitive to A^^^_|_^ and 

The four largest asymmetries for Oi are presented in fig. (3) as a function of cos(6'), 

employing the same cuts as used in 77 — > 77. At L — 500 fb~^ and Ai = 2 TeV, the 
"symmetric" asymmetries are statistically significant throughout the entire angular region, 
while the "antisymmetric" asymmetries are significant in the outer regions. In fig. (4) we 
present the total cross section versus ^/s for each inital polarization state, and the integrated 
asymmetry versus ^/s for the four symmetric initial states. Although a deviation in the cross 
section is not seen in any polarization state until high energies, the integrated asymmetry 
becomes quite large at low energies, suggesting that it provides a sensitive test of the anoma- 
lous couplings under consideration. So far we have only considered the operator Oj^; the 
distributions for the other two operators differ only in relative sign and overall magnitude 
at high energies. The asymmetries for polarization states (H 1 — ) and ( — h H — ) for all 

three operators are shown in fig. (5); those for OJ^ and are of similar magnitude but 
opposite sign. This can be understood in a simple way. At high energies, a transversely 
polarized Z effectively looks like a photon. When computing the amplitudes for the various 
operators in eq. (37), we must include all possible permutations of the three photons, as they 
are identical particles (for example, denoting the field strength tensors of the three photons 
by 1,2, and 3, the six permutations required to compute Oj^ become {12){3Z), (21)(3Z), 
(13)(2Z), (31)(2Z), (23)(1Z), and (32)(1Z)). The permutations required for Of and Of 
are identical to the permutations necessary to compute the amplitudes for the four photon 
operator of eq. (28); we can verify using the explicit forms of the amplitudes in the Appendix 

that the sum of Mj^iXi, A2; A3, A4) and M.2^{Xi, A2; A3, A4) gives ^^^^(Ai, A2; A3, A4), up to 
an overall constant. Since T violation leads to vanishing of the amplitudes M.±±±±, -^i^iqr 
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Table 1: Coefficients relating the C>ft^(2)xc/(i) to the Of. 



and A1±^=p± in 77 — > 77, the asymmetries for Oj and O2 must therefore be of opposite 
sign. SU(2) X U(l) operators whose coefficients 6^, 6^ are of the same sign will see destruc- 
tive interference between Of and Of, while those whose coefficients have opposite signs 
will see constructive interference. We will see similar behavior when investigating 77 ZZ. 

Table (1) presents the coefficients 6^ relating the C'^t^(2)xc/(i) ^ ^^e Of. The 7Z 
operators will interfere destructively in 0(^bb)(bb), 0(^ww)iww), a-nd 0(^BW)iBW) and con- 
structively in 0(^BB)(ww) and To estimate the value of A^, that can be probed 
at a 77 collider we have performed a combined least-squares fit to the normalized binned 
cross section, binned asymmetry, and total cross section, with ^/s — 1000 GeV for the two 

polarization states (-1 1 — ) and ( — h -I — ). These two choices are chosen to illustrate 

the sensitivities obtainable from both symmetric and asymmetric initial polarizations; the 
search reaches from the remaining polarization states are similar. The results are presented 
in fig. (6). Here, we have included only five of the SU(2) x U(l) operators; 0(^bb){ww) and 
0(ww){BB) differ only in the sign of their asymmetries, as do Obwbw and Owbwb, and 
hence have identical discovery regions. Remembering that our operators are of dimension 
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eight, and that the anomalous amphtudes therefore scale as ~ s^/A^, we see that the 
process 77 —>■ 'jZ is quite sensitive to the CP violating operators whose asymmetries inter- 
fere constructively. The attractiveness of this process is enhanced by the ease with which it 
can be experimentally reconstructed. 

In summary, we find that this process is a sensitive probe of CP violation in the gauge 
boson sector. The asymmetry considered here would be useful for any anomalous amplitudes 
A^^l'izr and M±^:t:± which contain imaginary parts. 



6 77 ^ ZZ 

In this section we will discuss the measurement of CP violation in 7(/ci) +7(/c2) —* Z{pi) + 



Z{p2). A detailed description of the SM process is given in |jT3 



We first briefly review the constraints imposed upon the ZZ amplitudes by Bose 
symmetry and parity violation. Let Aj = 0, ±1 denote the Z polarizations. Bose symmetry 
implies 

M'^,l,^is,t,u) = i-l)'^-'^MZ,,^{s,u,t) = MZ,xA'',^^,t) , (39) 
while oddness under P implies 

MZ,,,{s,t,u) = {-l)'^-'^-'M^_l,_,^_,^{s,t,u) . (40) 

These relations lead to the vanishing of the amplitudes Ai^^^-^ and A4^^^^. However, 
unlike in 77 — 77, the amplitude M±±±± is non-zero. Explicit formulae for these amplitudes 
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Figure 3: "Primarily symmetric" asymmetries (top) and "primarily antisymmetric" asym- 
metries (bottom) for Of, with Ai = 2 TeV, L = 500fb-\ and y/s = 1000 GeV. The bars 
indicate the corresponding statistical error. 
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Figure 4: Total (polarized) counting rate (top) and integrated asymmetry (bottom) versus 
y/s for Oj^, with Ai = 2 TeV and L = 500 fb~^. The bars indicate the corresponding 
statistical error, as well as a 1 % luminosity uncertainty in the event rate. The solid curves 
represent the SM event rates. 
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Figure 5: Comparison of A^z for the polarization states (H 1 — ) (top) and ( — h H — ) 

(bottom), with Aj = 2 TeV and L = 500 fb~^. Here, the statistical error bars have been 
suppressed. 
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Figure 6: Search reach for bC^ for polarization state (H 1 — ) (top), and (— + H — ) (bottom). 

All of the quoted sensitivities are at the 95 % CL. 
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are given in the Appendix. 

The SM amphtudes for this process are similar to those for 77 — > 77 and 77 — > 7Z 
in that they are completely dominated at high energies by the imaginary parts of ■M±±±±, 
and Ai^^^^. However, they differ in one important aspect. At one loop they receive 
contributions from the Higgs sector through the one loop h'-ff vertex, and are therefore not 
independent of the mechanism of electroweak symmetry breaking. The contributions to this 
process from a neutral Higgs, such as the one present in the SM, have been given in [0]; 
for a Higgs mass nih < 2m z (which includes the value rrih = 115 GeV suggested by LEP), 
this contribution is negligible. For larger values of the Higgs mass, the two Z bosons will 
reconstruct to rrih for this contribution, making it separable from continuum ZZ production. 
We thus neglect the Higgs contribution in our analysis. 

The structure of both the SM and the anomalous amplitudes are similar to those in 
77 7Z; in particular, the relevant asymmetry that contains interference between the large 
imaginary SM amplitudes and the anomalous amplitudes is again given by 

[ da_\ _ ( da_\ 
dO.) \dQ. 

Varying the intial polarization state (Pel, -P/i, Pe2i P12) yields less information than in 77 — > 
7Z, however, for reasons discussed below. The numerator of the asymmetry again takes the 
form 

Azz oc {Uxi)+Ux2))[\M++\^ -\M- 

+ iU^i) - Ux2)) [\ -M+- P - I -M-+ \'\ , (42) 

but the vanishing of Ai^^^z^ and Ai'^^^^ implies that the second term in eq. (42) equals 
zero. We will only observe the "symmetric" asymmetries, and the "antisymmetric" asymme- 
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tries for the initial states ( — h H — ) and (+ H ) will vanish. We can again simultaneously 

measure the polarized differential cross section, 



da 



pol 



1 



dcos{9) 



2567rs 



l + Uxi)Ux2)){\M++ \' + \M- 



(43) 



The anomalous structures relevant for the interaction 77 ZZ can be found in eq. 
(20). We will write them in the now familiar form 



Of 



zz 



(A 



or 



{F^,F'''')iZ^^ Z 



pa J 
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zz 



(Aa)^ 



{F^,Z^'''){FP^ Z 



pa) 



o, 



zz 



O 



zz 



^ Z^'^F.pZP'' Fa, 



(A5) 



(44) 



We have again set gejf = 1, where geff is the coupling constant associated with each opera- 
tor. The detailed amplitudes are presented in the Appendix. As in 77 — 7Z, those with one 
longitudinal Z are suppressed by a factor of rriz / y/s, and we will not attempt to reconstruct 
them. Amplitudes with two longitudinal Zs are 0(m|/s), and hence negligible in our ap- 



proximation, which again entails the use of the asymptotic SM expressions presented in [0]. 
The amplitudes M±±±± vanish for , and therefore this operator does not contribute to 
Azz- 
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The experimental reconstruction of this process is more difficult than for 77 — > 77 
or 77 — >■ 7Z, as there is no high photon to tag. Event reconstruction will likely require 
demanding a leptonic Z decay or possibly a Z ^ bb decay for one of the Z bosons, resulting 
in a detection efficiency of roughly 40 ~ 45%. However, it is beyond the scope of this paper 
to present a detailed ZZ reconstruction analysis; such a study should also make use of the 
complete SM amplitudes and detector efficiencies. An appropriate scaling of the integrated 
luminosity when observing the quoted sensitivities is sufficient to estimate the effects of 
imperfect reconstruction. 

The binned asymmetries for the symmetric initial polarizations of Of^ are shown in 
fig. (7). They are typically larger than those for 77 —>■ jZ. Note that the central value 
of the (H ) asymmetry displays a peculiar peaked shape. The integrated asymme- 
tries and total cross sections for the same initial polarization states are shown in fig. (8). 
Again, the integrated asymmetry becomes large at low energies for all initial polarization 
states, suggesting that this asymmetry is a sensitive test of these anomalous operators. The 
asymmetries for Cf'^ are of similar magnitude as those for Of^ , but have the opposite 
sign. Similarly, the asymmetries for Of^ and O^^ possess opposite signs, and are somewhat 
smaller in magnitude than those for Of^. SU(2) x U(l) operators whose coefficients c^, 
have the same sign will exhibit destructive interference between Cf^ and Cf ^, while those 
with coefficients of opposite signs will see constructive interference; the same interference 
pattern exists for the coefficients and c^. 

To convert these results into statements regarding the sensitivity of this process to 
the various SU(2) x U(l) operators we need the coefficients presented in table (2). Unlike 
the situation in 77 — > 7Z, we see that all of the interference effects are destructive, given 
the relative signs of the for each operator. However, since c^r ~ 0.59 and ~ 0.053, a 
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Table 2: Coefficients ^ relating the Cff^(2)xc/(i) to the 0| 



complete cancellation does not occur except for the three operators 0(^bb){bb), 0(^[yw){ww)-, 
and 0(^BW){BW)- We have performed a combined least-squares fit to the normalized binned 
cross section, binned asymmetry, and total cross section, with ^/s = 1000 GeV, for the 

polarization state (H 1 — ) to estimate the value of Aq, that can be probed in this process 

at a photon collider. The search reaches obtainable from the other symmetric polarization 
states arc similar. The results arc presented in fig. (9). We only display results for the 
operators 0^bb){bb), 0(bw){bw), 0(bb){ww), and Obwbw', the sensitivity to 0(bb){bb) and 
0{ww){ww) is identical at high energies. Similarly, 0(^bb){ww) and 0{ww){bb) differ only in 
the sign of their asymmetries, as do Obwbw and Owbwb-i and hence yield identical search 
reaches. We see that overall, this process is not quite as sensitive to the anomalous operators 
as 77 ^ 7Z due to the destructive interference arising from the SU(2) x U(l) embedding 
of the 77ZZ vertex structures. However, those operators for which the asymmetry does not 
vanish can still be constrained quite stringently. 

In conclusion, the sensitivity of Azz to interference with the large SM amplitude 
Ai±±±± renders this process a sensitive probe of CP violation appearing in gauge boson self 
couplings. 
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7 Conclusion 



In this paper we have examined the possibihty of using the processes 77 77, 7Z, and 
ZZ to search for sources of CP violation that manifest themselves at energy scales above 
those of future colliders. We have parameterized these effects by an effective Lagrangian 
containing the complete and independent set of CP-odd, SU(2) x U(l) invariant quartic 
gauge boson operators relevant to these interactions. The considered processes vanish at 
tree level in the SM, and at one loop are completely dominated by the imaginary parts of 
the amplitudes A4±±±±, ■M±^±^, and M.±z^^±. The fact that they might be powerful tools 



in indirect searches for new physics was pointed out in [11], |T2|, |T3[ , where their sensitivity 
to supersymmetric loop effects was discussed. Our results show that they are also of utility 
in searches for sources of CP violation that contribute to self interactions of gauge bosons. 

To demonstrate that these processes are indeed sensitive probes of the gauge boson 
sector, we compare our results with four others found in the literature. We first note that 
while our operators are of dimension eight, those with which we compare are of dimension six. 
Although our comparisons are for the search reaches for the new physics scale A, the reader 
should be aware that a given limit on A translates into a greater sensitivity to a dimension 
eight operator, which scales as s^/A^, than a dimension six operator, which scales as s/A^, 
for ^/s < A. We also note that the constraints on the anomalous SU(2) x U(l) operators 
obtained in this paper are stronger than those implied by unitarity. We first contrast our 
search reaches with the direct bounds on quartic gauge boson couplings given recently by the 
LEP experiments [jl8|. They consider operators of the form (e^Oj/lGA^) O, and obtain bounds 
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Figure 7: Binned asymmetries for symmetric initial polarizations for Of ^, with Ai = 2 TeV, 
L = 500 fb~^, and y/s = 1000 GeV. The bars indicate the corresponding statistical error. 
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Figure 8: Total (polarized) counting rate (top) and integrated asymmetry (bottom) versus 
^/s for Of^, with Ai = 2 TeV and L = 500 fb~^. The bars indicate the corresponding 
statistical error, as well as a 1 % luminosity uncertainty in the event rate. The solid curves 
represent the SM event rates. 
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Figure 9: Sensitivity to for the polarization state (H \ — ), at the 95 % CL. 

of a,/16A2 < IQ-^GeV^ on CP-conserving iy+iy~77, W+W~Z-f, and vertices. The 

operators considered in this paper are of the form (e^/A'*) O (with Qeff = 1), and therefore 
our search reach for 1/A^ are to be contrasted with their bounds on aj/16A^. As expected, 
our sensitivity at a higher energy and higher luminosity collider are very significantly better. 
We next compare our results with those of a study of these same operators at a 500 GeV 
e+e~ collider with an integrated luminosity of 500 fb^^ ||19|. In this study the operators are 
parameterized in the form (Ke'')/(4s^c^A^) C, and bounds of (Ke^)/(4s^c^A^) < 0.13 x 

10"^ GeV^ obtained on vertices such as the ZZ77 vertex. Setting k = 1, this results in 
A > 800 ~ 900 GeV, at the 3a level. The sensitivities obtained in this paper, normalized to 
the appropriate luminosity and energy are A > 800 ~ 1700 GeV, where the better constraints 
are for those operators contributing to one of the considered asymmetries. The search reaches 
obtained here on A, after adjusting the statistical significances, are still stronger, even though 
the operators we consider are of dimension eight. The sensitivity of future hadron colliders 



such as the LHC to these anomalous CP-even 77ZZ and 'yyW^W has been studied in [^] 
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This paper obtains search reaches of A > 800 ~ 1100 at the 95 % CL, which are again 
somewhat weaker than those we have presented. Finally, the authors of P] consider CP-odd 
dimension six operators contributing to the process 77 —>■ W~^W~ at ^/s = 1 TeV assuming 
20fb~^ of integrated luminosity. They cast their operators in the form {Y/ 171^)0, and 
obtain bounds of F < 2 x 10^'^ ~ 4 x 10^^. Note the we have improved their results by a 
factor of 2.5 to simulate the search reach assuming 50fb~^ of integrated luminosity. These 
translate into limits of A > 500 ~ 4100 GeV in our notation, at the la level. The relevant 
limits in this paper to compare to are those assuming 50fb~^ of integrated luminosity, which 
range from A > 800 GeV to A > 2500 GeV at the 95 % CL. Our constraints on the new 
physics scale A are not quite as strong, which is not too surprising considering the large cross 
section for 77 — > W^W~ and the fact that their operators are of dimension six while ours 
are of dimension eight, although the differences are lessened after accounting for the greater 
statistical significance of our results. However, the processes considered are complementary; 
they study the charged gauge boson production process 77 — W^W~ while we examine the 
neutral gauge boson processes 77 77, 7Z, and ZZ. 

In summary, the processes 77 — >■ 77, 7^, and ZZ are sensitive probes of CP violation 
in the gauge boson sector. The sensitivities obtained here to the anomalous CP-odd oper- 
ators affecting these interactions compare favorably with similar limits on anomalous gauge 
boson couplings found in the literature. In addition, the examination of these processes 
nicely complements previous studies that have focused primarily on W boson, top quark, or 
Higgs production. Their utility in other indirect physics searches at future photon colliders 
should certainly be studied in more detail. 
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Appendix 

We collect here formulae describing the various photon distribution functions; deriva- 
tions are given in M. We first define the auxiliary functions 



1 



1 — X 



+ (1 - x) - 4r(l -r)- P^Pir z{2r - 1)(2 - x) 



(45) 



where r = x/[z{l — x)], and 
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zj ^ ^ 2 z+1 2{z + lf 



(46) 



Here 2; is a variable describing the laser photon energy, and is given by 2; = AEeEi/ml, 
where E^. is the electron beam energy and Ei the laser energy. Varying z changes the 
maximum value of the backscattered photon beam energy Xmax, where Xmax = z/{l + z). 
We will set z = 2(1 + V2), which maximizes Xmax while preventing interactions between the 
backscattered photons and laser beam. In terms of these functions and variables the photon 
number and helicity distributions take the form 
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/ is the photon number density function, while and are respectively the circular and 
linear helicity distribution functions. We can now write the observable differential cross 
section as 

where s = AE^ and n\ accounts for any final state Bose symmetries. 

Here we present the kinematics and amplitudes used for the processes in this paper. 
The generic interaction considered is Fi{ki, ei) + -^2(^:2, €2) — > Vi{pi, Vi) + V2(p2, V2), where 
Fi, F2 are the incoming photons and Vi, V2 the outgoing gauge bosons. The momenta and 
polarization vectors take the following form in the c.m.s. frame: 

= ^(1,0,0,1) k2i, = ^(1,0,0,-1) 

pi^ = ^il,Pse,0,Pce) P2^ = ^(l,-/3se,0,-/3ce) 



^1 = -;j5(0,Ai,^,0) e^2 = ^=,{0, -\2,t,0) (49) 

V^L = ^{-(3,se,0,ce) V^l = -^(/3, se, 0, c,) 

with sq — sin^, cq — cos9, Pv — y^l — 4my/s, and s — {ki + k2Y — (pi +^2)^- Given these 
expressions, it is straightforward to compute the amplitudes for the operators of eq. (20). 
As mentioned in the text we neglect terms of 0(m^/s, m|/s), as the SM amplitudes we use 
also omit terms of this order. We find that this approximation is numerically valid for the 
energies considered here. The amplitude for the anomalous 77 operator of eq. (20) takes 
the form 

MViK A2; A3, A4) = [(A3 + A4)(l + A1A2) - (Ai + A2)(l + A3A4)J (3 + c^) . (50) 

This amplitude is exact, as mz and mw appear only in the SM amplitudes. Note that 
this amplitude is odd under the interchange (Ai, A2) (A3, A4); this is a consequence of T 

36 



violation and the fact that the two final state particles are the same as the initial states. 
The anomalous amplitudes for the process 77 — > 7Z (with V2 = Z) are, neglecting terms of 
Oiml/s), 
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(51) 



Again neglecting terms of 0{m\/s), the ZZ amplitudes become 



A^r(Ai,A2;A3,A4) = -^^^pyi (A3 + A4) (1 + A1A2) 
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32(AP)4 



(1 + A1A2) - A3(Ai + A2) 
6(Ai + A2)(l + A3A4) + 2(A3 + A4)(l + A1A2) 



+2 - (Ai + A2)(l + A3A4) + (A3 + A4)(l + A1A2) 



lerazs s 
4(AP)4 V 2 



(I + A1A2) - A3(Ai + A2)ce 
6(A3 + A4)(l + A1A2) - 2(Ai + A2)(l + A3A4) 



+2UA3 + A4)(l + A1A2) - (Ai + A2)(l + A3A4) jc^ 



lemzs Is 
4(AP)4^' V 2 



(I + A1A2) - A3(Ai + A2)ce 



(52) 



Amplitudes with two longitudinal Z bosons are of 0(m|/s), and hence are negligible in our 
approximation. 
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